Dual fermion approach to nonlocal correlations in the Hubbard model 
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A new diagrammatic technique is developed to describe nonlocal effects (e.g., pseudogap forma- 
tion) in the Hubbard-like models. In contrast to cluster approaches, this method utilizes an exact 
transition to the dual set of variables, and it therefore becomes possible to treat the irreducible 
vertices of an effective single-impurity problem as small parameters. This provides a very efficient 
interpolation between weak-coupling (band) and atomic limits. The antiferromagnetic pseudogap 
formation in the Hubbard model is correctly reproduced by just the lowest-order diagrams. 
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Lattice fermion models with a strong local interaction 
(Hubbard-like models [1]) are believed to catch the basic 
physics of various systems, such as high- temperature su- 
perconductors |2|, Isj , itinerant-electron magnets [J, Mott 
insulators [s^], ultracold atoms in optical lattices [y], etc. 
Unfortunately, the analytical treatment of these prob- 
lems is essentially restricted by the lack of explicit small 
parameters for the most physically interesting interac- 
tions. Direct numerical methods, such as exact diago- 
nalization p"] or quantum Monte Carlo (QMC) |8|, l9| are 
limited by the clusters being of a relatively small size, 
or face other obstacles such as the famous sign problem 
for QMC simulations at low temperature [lO[. There is 
a very successful approximate way to treat these models 
via the framework of so-called dynamical mean-field the- 
ory (DMFT) [5], where the lattice many-body problem 
is replaced with an effective impurity model. This ap- 
proach is essentially based on the assumption of a local 
(i.e. momentum-independent) fermionic self energy. In- 
deed, there are numerous interesting phenomena which 
are basically determined by local electron correlations, 
such as Kondo effect (uj, Mott-Hubbard transitions Q 
and local moment formation in itinerant-electron mag- 



nets [12|. At the same time, momentum dependence of 
the self energy is of crucial importance for Luttinger liq- 
uid formation in low-dimensional systems [j, Il3i] , d-wave 
pairing in high- Tc superconductors [1, Q, EH , and non- 
Fermi-liquid behavior due to van Hove singularities in 
two dimensions [3]. Recently a rather strong momen- 
tum dependence of the effective mass renormalization in 
photoemission spectra of iron was observed pj] . 

Currently, non-local many body effects in strongly cor- 
related systems are mainly studied via the framework of 
various cluster generalizations of DMFT 13, Hi, [1^. 
Cluster methods do catch basic physics of d-wave pair- 
ing and antiferromagnetism in high- Tc superconductors 
14, lij, and the effects of intercite Coulomb interaction 
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in various transition-metal oxides [2C 
same time, however effects like Luttinger liquid formation 
or van Hove singularities can not be described in cluster 
approaches. In such cases the correlations are essentially 



long-ranged and it is more natural to describe the cor- 
relations in momentum space. Recently attempts have 
been made to consider non local correlation effects in 
momentum space starting from DMFT as a zeroth-order 
approximation 2^, l^l • This approach requires a solution 



of ladder-like integral equation for complete vertex T and 
the subsequent use of the Bethe-Salpeter equation to ob- 
tain Green's functions. The first step here exploits an 
irreducible vertex of the effective impurity problem 7^^\ 
whereas the second step uses just the bare interaction 
parameter U . This second step makes the generalized 
DMFT approoach mostly suitable to the weak-coupling 
regime [25]. 

In this Letter, we present a scheme which is accurate 
in both small- /7 and large- limits and does not require 
numerically expensive solutions of any integral equations. 
A comparison of the results with lattice QMC simulations 
for the two-dimensional (2D) Hubbard model in the pseu- 
dogap regime demonstrates that the scheme is actually 
accurate even in the less-favorable case of intermediate 
U. 

We proceed with 2D Hubbard model with the corre- 
sponding imaginary-time action 



S[c,c*] 



E 

(joka 



(1) 

Here /3 and /i are the inverse temperature and chemical 
potential, respectively, co = {2j + l)7r//3, j = 0, ±1, ... are 
the Matsubara frequencies, a =t, | is the spin projection. 
The bare dispersion law is Ck = — 2t(cos kx-\-cos ky)^ c* 
are the Grassmannian variables, n^crr = c*^r^«crr, where 
the indices i and k label sites and quasi-momenta. 

In the spirit of DMFT, we introduce a single-site refer- 
ence system (an effective impurity model) with the action 

^imp = y](Aa; -/i-zcj)c*^^c,^,^ + [/ / n^rnirdr (2) 

where A^^ is as an yet undefined hybridization function 
describing the interaction of the effective impurity with 
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a bath. We suppose that ah properties of the impurity 
problem are known, so that its single-particle Green's 
function is known, and the irreducible vertex parts 
7^^\7^, etc. Our goal is to express the Green's function 
Gouk and vertices T of the lattice problem in Eq.([T]) via 
these quantities. 



Since A is independent of /c, the action ([T]) can be 
represented in the form 

= ^6'imp[ci,c*] - ^(A^ - ek)cl^^Cujka- (3) 

i ioka 
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We utilize a dual transformation to the set of new Grassmannian variables /, /* . The following identity 

is valid for arbitrary complex numbers A and B. We chose = (A^^ — Ck) and B'^ = g~'^{Au; — e^)"^ for each set of 
indices uj^k^a. 

With this identity, the partition function of the lattice problem Z = J e~^^^'^ ^Vc*Vc can be presented in a form 
Z = Zf j j e-^[^'^*"^"^*]p/*P/Pc*Pc, where 

i ooka 



and Zf is a product U^kgli^uj - ^k)- 

As a next step, we establish an exact relation between 
the Green's function of the initial system Gr-r'.i-i' = 

- < TcriC%>^, > and that of the dual system GfTS,i-i^ = 

— < Tfrif^f^f >. To this aim, we can replace Ck 

^k + ^^uk with a differentiation of the partition function 
with respect to Se^k- Since we have two expressions for 
the action ([1]) and ([5]), one obtains 



^(A^-6,)-^G^T^ + (A^-6,)-\ (6) 



where the last term follows from the differentiation of Zf. 

The crucial point is that the integration over the ini- 
tial variables c* , ci can be performed separately for each 
lattice site, since Y.k ifk^k + c*/^) = ifi^i + 
For a given site i, one should integrate out c*,q from 
the action that equals 5'si^e[ci, <, /i, /*] = 5'imp[Q, <] + 
Eo. ^j'(/:cc. + c^/c). We obtain 



(7) 

where Zimp is a partition function of the impurity prob- 
lem (|2]). The Taylor series for V[fi, /*] in powers of /i, /* 
starts from —71234/1/2/3/4 (indices stand for a combina- 
tion of (J and cj, for example fl means f^^^^^^)- Further 
Taylor series terms yeld similar combinations including 
^(^) of higher orders. 

We arrive with an action S depending on the new vari- 
ables /, /* only: 



s[f. /*] = E ((^- - ^^)"' + 3.) r.kaUka + E 

i 

(8) 
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with Vi = V[f*^fi]. In this dual action, the interac- 
tion terms remain localized in space, but are non-local in 
imaginary time, since, for example 7^^^ depends on the 
three independent Matsubara frequencies. To obtain the 
dual potential V for a practical calculation, one should 
solve then the impurity problem ([2]). 

Finally a regular diagrammatic expansion in powers 
of V can be performed. We draw skeleton diagrams, so 
that the lines in diagrams are renormalized dual Green's 
function, whereas the vertices are 7^"^^. The rules of di- 
agram construction are very similar to usual ones, but 
the six-leg and higher-order vertices appear because 7*^^^ 
and higher terms are present in the series for V. Figure 
1 shows several diagrams contributing dual self-energy 

We use the skeleton-diagram expansion for the dual 
self-energy since it leads to the conserving theories, ex- 
actly like in conventional diagram technique 27, 2^. 
The Baym criterion of a conservative theory is the ex- 
istence of a functional of the Green function ^[G] such 
that 1^ = S. Here, the variation SG comes from the 
infinitesimal variation of the Gaussian part of the ac- 
tion, 6{G^)~^c*c. In our consideration, we consider 
also the infinitesimal variations of the dual potential 
^i^duai)~^^*^' ^^^^ approximation dually 

derivable, if there exists a functional ^duai^Qduai^ such 

as ^Qduai = where the variation comes from 

^(^duai)~^ ' Now, it turns out that the theory is ^- 
derivable if it is dually <l>- derivable. The proof uses the 
relation between functional <l> and the partition function 
In Z = - TrSG ;^r In(-G) + C (here G is an additive 
constant; see Ref.27, Eq.(47)). Since a similar relation 
takes place for and InZ and since the partition 
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FIG. 1: (color online) Various diagrams for T^"^^"^^ and the 
scheme of calculation. The calculation includes "big" and 
"small" loop, marked with red and black lines, respectively. 
The small loop is to determine the renormalized dual Green's 
function G'^^^^ in a self-consistent way, for given A,^, and 
^(^) ]3jg Jqqp jg determine A. Only the big loop 

requires a solution of the impurity problem. 

function is the same for the initial and dual variables, 
this gives a one-to-one correspondence between ^[G] and 
^duai^Qduai^ xhis is just a sketch; the detailed proof will 
be published elsewhere. 

It is important to understand what can be a small pa- 
rameter in the expansion in dual diagrams. Clearly, if U 
is small, then 7*^^^ ex [/, 7*^^^ ex U'^ etc., and in the weak- 
correlated regime vertices in the diagrams will be small 
(Fig. 1), and higher-order vertices will be even smaller. 

At this point we establish a condition for A, which 
was so far an arbitrary quantity. We use a self-consistent 
condition 

k 

It means that the simple closed loops in diagrams van- 
ish. In particular, this leads to the vanishing of the first- 
order "Hartree" corrections in the diagrammatic expan- 
sion. The diagram series of this kind has several im- 
portant peculiarities. First of all, let us consider the 
zeroth-order approximation, Y^^^^^^) = 0. In this case, 
the condition (|9]) becomes 

y rT=0- (10) 



It is easy to show that this is exactly equivalent to the 
DMFT equation for the "hybridization function" A^^ ^ . 
It is known that DMFT behaves correctly near the atomic 
limit. In terms of the dual variables, one can observe 
that since e, A <C near the atomic limit, it follows 
from the condition ([10|) that G^^^^ ^ glck « guj in this 
case. It easy to check that this argumentation is valid the 
scheme of an arbitrary diagrammatic order: the disper- 
sion of G^^^^ is small near the atomic limit and therefore 
([9]) means that lines in dual diagrams carry a small factor 
eg~^. This ensures the fast convergence of new diagram- 
matic expansion in the strong-coupling limit. 

As the most challenging test of our approximation 
scheme in the intermediate regime, we performed the cal- 
culation for the half-filled square-lattice Hubbard model, 
at sufficiently low temperature /3~'^ = \t\/5. The value 
of U was varied from small numbers to a bandwidth 8|t|. 
The block scheme of our calculation is shown in Fig. 1. 
It has a good practical convergence: typically, about 10 
iteration are enough to ensure convergence. 

In order to obtain reference point for a further compar- 
ison with the results of our new approximation scheme, 
we performed a direct lattice QMC calculation with the 
continuous-time QMC code [30|. There are strong anti- 
ferromagnetic fluctuations in the system, although true 
antiferromagnetism is impossible at finite temperature 
in the 2D system with an isotropic order parameter [2^ . 
Consequently, the increase of U results in a formation of 
an antiferromagnetic pseudogap. 

It was also noticed that single-site DMFT calculation 
for this system shows no pseudogap in the density of 
states, although the data for local part of self-energy are 
reproduced quite well in DMFT. Thus we concluded that 
the formation of pseudogap is entirely related to the non- 
local part of S, neglected in DMFT. 

We present the results of the dual-fermion calculation 
with only diagram (b) taken into account. All other dia- 
grams are smaller both in the strong-coupling and weak- 
coupling regime, due to extra vertices or extra lines, re- 
spectively. Computational results are illustrated by Fig. 
2. The upper panel shows an imaginary part of the 
self-energy. In the DMFT this quantity is momentum- 
independent. Our calculations show a very strong k- 
dependence with a maximum near the Fermi surface. At 
relatively small value U = 1 the peaks of ImS are located 
near the van Hove singularities (left picture), as it can be 
understood from the weak-coupling expansion. Contrary, 
for an antiferromagnetic system near the atomic limit, 
lmTik,uj=o would be a simple delta-function peaking at 
Fermi surface. For a pseudogap regime at finite /7, the 
width of this peak is of course finite, but the altitude al- 
most does not depend on the point at Fermi surface (right 
picture). The lower panel shows an effective renormalized 
dispersion law Ck + ReT^k,uj=0' For the metallic regime, 
the renormalization is small. For an antiferromagnetic 
insulator, the would be a pole in IleT^k,u^o at the Fermi 
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FIG. 2: (color online) Fermi-energy properties of the half- 
filled Hubbard model calculated with the leading dual dia- 
gram correction b. The calculations have been performed for 
the bandwidth 8t = 2 at /3 = 20, for different values of U. 
Upper panels are contour plots for ImS/c at Fermi energy. 
At U — 1, ImS peaks in the four van Hove points, whereas 
lm^{U = 1) is approximately constant in all points of the 
Fermi surface. Note also that the change from [/ = 1 to 
U — 2 leads to a 10^ increase in ImS. Lower panel shows 
a graph of the effective dispersion law, Ck + ImS/c at Fermi 
level, plotted along the (0,0) — (7r,7r) direction. The ini- 
tial "cosine" dispersion law Ck is almost not renormalized at 
[7=1. Contrary, for U — 2 the curve shows the antiferro- 
magnetic properties. The result of direct QMC 10 x 10 lattice 
simulation are shown with dots and confirm this picture. 



surface. For the pseudogap regime, fluctuations virtually 
move this pole from the real-frequency axes, as the curve 
for U = 2 shows. 

Thus, our scheme continuously interpolate between the 
two very different regimes. It should be stressed that the 
quantities under study have very strong /c-dependence 
and that it would be very difficult to obtain the result of 
this kind, for example, in cluster calculations. Whereas 
for the weak-coupling regime effective schemes to calcu- 
late nonlocal self-energy are known, such as FLEX [2^ 
or parquet to our knowledge, there is no alternative 
scheme yet for the strong coupling case. 

To conclude, we have formulated an effective pertur- 
bation theory to calculate the momentum dependence of 
self energy starting with single-site DMFT or any local 
approximations. The vertices of the effective impurity 



problem play the role of formal small parameters. Due 
to the transformation to dual fermionic variables, consid- 
eration of a few leading diagrams provides a quite satis- 
factory description of the nonlocal correlation effects in a 
broad range of parameters, up to the atomic limit. This 
scheme can be easily generalized to multiband case to 
be implemented into realistic electronic structure calcu- 
lations for strongly correlated systems. 
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